Microscopic properties of vortex state in YBe probed by muon spin rotation 
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Local magnetic field distribution B(r) in the mixed state of a boride superconductor, YBe, is 
studied by muon spin rotation (/iSR). A comparative analysis using the modified London model 
and Ginzburg-Landau (GL) model indicates that the GL model exhibits better agreement with /xSR 
data at higher fields, thereby demonstrating the importance of reproducing the field profile near 
the vortex cores when the intervortex distance becomes closer to the GL coherence length. The 
temperature and field dependence of magnetic penetration depth (A) does not show any hint of 
nonlocal effect nor of low- lying quasiparticle excitation. This suggests that the strong coupling of 
electrons to the rattling motion of Y ions in the boron cage suggested by bulk measurements gives 
rise to a conventional superconductivity with isotropic s-wave pairing. Taking account of the present 
result, a review is provided for probing the anisotropy of superconducting order parameters by the 
slope of A against field (77 oc dX/dH). 



I. INTRODUCTION 

The revelation of superconductivity in magnesium di- 
boride (MgB2, Tc — 39 K) has stimulated renewed in- 
terest to other boride superconductors as an alterna- 
tive path to novel superconductors with ever higher Tc.^ 
Among those, yttrium hexaboride (YBe) currently holds 
the position of second highest transition temperature 
{Tc > 7 K),^ thereby drawing considerable attention on 
the detailed mechanism of superconductivity. In partic- 
ular, it stands as a focus of investigation on the effect of 
rattling motion of alkaline earth metals contained in a 
relatively large boron cage that gives rise to anharmonic 
phonons having relatively low energy (c^ 10^ meV).^ A 
similar situation is reported in a variety of novel super- 
conductors including metal-doped silicon clathrates^ and 
/3-pyrochlores.^ 

While YBe is known as a type II superconductor in 
which the Ginzburg-Landau (GL) parameter (= A/^, 
with A being the magnetic penetration depth and ^ = ^gl 
the GL coherence length) is greater than unity, it is char- 
acterized by a low upper critical field [Bc2 = /io^c2(0) ^ 
0.3 T] and associated long ^ (~ ^0 — 30 nm, the BCS 
coherence length) that is only a few times smaller than 
A (^ 130 nm). This means that the volume fraction 
occupied by vortex cores {c:^ 7r^^i^/<I>o, where H is the 
external magnetic field and <I>o is the flux quantum) is rel- 
atively large in the flux line lattice (FLL) state over the 
entire field range, making YBe an ideal stage for study- 
ing the electronic structure of vortex cores in detail. In 
this situation, it is pointed out that the modified London 
(m-London) model, which is known to be an excellent 
analytical model to describe the spatial distribution of 
magnetic field B{r) for arbitrary Hi but at low magnetic 
induction {B/Bc2 < 0.25),^^^ may not be appropriate for 
the basis to extract physical parameters out of B{t).^ 



One of the primary purposes of the present study is to 
address this potential problem by making a comparative 
analysis of /iSR data by both m-London and GL mod- 
els, which reveals that the latter provides systematically 
better description of magnetic field profile at higher mag- 
netic induction. Thus, it is inferred from this study that 
we have to resort to the GL model for proper understand- 
ing of the FLL state over the wider field range. Fortu- 
nately, however, in contrast to the case of classical type 
II superconductors such as Nb or V (a^ :^ 1) where one 
might need numerical approach to obtain B{r) with suffi- 
cient precision^ii^ the present result suggests that an an- 
alytical solution of the GL model obtained by variational 
methodiii^ is sufficient for the system with reasonably 
large hi. 

It must be stressed that the magnetic penetration 
depth controlling the actual shape of B{t) corresponds 
to an effective value of the London penetration depth 
determined by the superfluid density n^. 



(where m* is the effective mass of the charge carriers), 
so that it may vary according to the change in Ug due to 
quasiparticle excitations (breaking of the Cooper pairs) 
of various origins. In particular, it is anticipated that the 
pair breaking due to the quasiclassical Doppler shift of 
the Fermi momenta around vortices would lead to the low 
energy quasiparticle excitation (the nonlinear effect) 
Considering that the number of vortices is proportional 
to the external field, it is predicted that Ug is more re- 
duced at higher H due to the Doppler shift. Another 
important consequence of the anisotropic order param- 
eter is the nonlocal effect due to the variable coherence 
length [^0 = ^o(k)] that may exceed the London pen- 
etration depth (Al) over a certain region of the Fermi 
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surface, as it is inversely proportional to the order pa- 
rameter [^o(k) cx fi'UF/7rA(k), with vp being the Fermi 
velocity]. In this situation, the flow of supercurrent is 
strongly modified over the region Al < ^o(k), leading to 
the change (expansion) of effective A.^^ Note that this is 
not limited to the case of nodal gap, but would be in ef- 
fect for anisotropic gap (or multi-gapped) with the mini- 
mum gap satisfying Amin < ^^f/ttAl- Meanwhile, multi- 
gapped superconductors have a multitude of coherence 
length and associated vortex core radius within which the 
quasiparticles are confined. The presence of quasiparti- 
cles over the region of larger vortex cores would serve as 
a factor to enhance the effective A more strongly at lower 
fields, thus mimicking the above mentioned Volovik effect 
without gap nodes. 

We have shown that the magnetic field dependence 
of A provides a useful criterion to assess the anisotropy 
of superconducting order parameter; the greater field 
gradient {dX/dH) corresponds to stronger anisotropy or 
multi-gapped structure in the order parameter We 
demonstrate in YBe that A deduced from the GL model- 
based analysis does not depend on the external magnetic 
field. This provides strong microscopic evidence for the 
isotropic superconducting gap in YBe. 



II. MODELS FOR THE FIELD PROFILE 

In the conventional /iSR studies of the FLL state, pos- 
itive muons are implanted to the superconducting speci- 
men with their initial spin polarization perpendicular to 
the external magnetic field [P(0) -L with J? || i in the 
following], so that they may probe B{t) [= ^^(r)] gener- 
ated by a perio dic array of magnetic vortices which are 
apart by ao ^o/H sinO (with the apex angle = 60° 
for the hexagonal FLL and = 90° for the square FLL). 
Considering that muon stops at the interstitial sites in 
the atomic unit cell of the crystalline lattice which has a 
periodicity much shorter than that of vortices (ao > 89 
nm for /hqH < 0.3 T assuming a hexagonal FLL), one 
can presume that /iSR signal in the FLL state provides 
a random sampling of 5(r), 

P{t) = P,{t)^iPy{t) 

/oo 
n{B) exp{i-f^Bt - icfydB, (2) 
-co 

n{B) = {S{B - B{v))),, (3) 

where n{B) is the spectral density for the internal field 
defined as a spatial average (()r) of the delta function, 7^ 
is the muon gyromagnetic ratio (= 27r x 135.53 MHz/T), 
and (/) is the initial phase of rotation.^ These equations in- 
dicate that the real amplitude of the Fourier transformed 
muon spin precession signal. Re / P{t) exp{—iji^Bt)ji^dt^ 
corresponds to the spectral density, n{B). 

In general, B{t) is obtained by solving the GL equation 
for the field profile and the order parameter ip{r) in the 
self-consistent manner, which usually requires numerical 



approach. However, in the case of strong type II super- 
conductors {k, ^ 1), the nonlocal character is confined 
within the immediate vicinity of vortex cores (|r| < ^gl) 
and the London approximation can be extended by intro- 
ducing a cutoff factor associated with the vanishing ?/^(r) 
at the vortex center. In the modified London model, B{r) 
satisfies the equation 

B{t) - curl[A • curl5(r)] = ^0 X] p(r - r^), (4) 



where A is the tensor related to the penetration depth 
and p(r) is the source term for the vortices located at 
the position in the lattice. The solution is given by a 
sum of the magnetic induction from isolated vortices to 
yield 



B{t) = ^6(K)exp(-zK-r) 



K 



6(K) 



Bo 



1 + KiK 



(5) 



(6) 



where K are the vortex reciprocal lattice vectors, Bq 
{0^ jiioH) is the average internal field, I/(K) is the effective 
London penetration depth, F{K^ ^c) is the cutoff factor 
(with (,c oc being the cutoff parameter) to remove the 
unphysical divergence in the magnetic field distribution 
near the vortex center: The local London approximation, 
p{r) = (^(r), corresponds to F{K^^c) = 1, where B{r) ex- 
hibits a logarithmic divergence. The local approximation 
also leads to the reduction of the term KLK into X^K'^ 
corresponding to the isotropic Fermi surface. Here, "lo- 
cal" implies that the electromagnetic interaction does not 
depend on the wave vector (K), whereas that referring to 
the source term is primarily meant for = 0. The form 
of F(i^, ^c) depends on the modelling of p(r) and asso- 
ciated ?/^(r) in the vicinity of vortex cores. For example, 
the Gaussian cutoff. 



F(i^,ec) =exp 
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(7) 



which is widely used for /iSR data analysis (valid for 
B/Bc2 < 0.25 and />: ^ 1), is derived by the isotropic 
GL theory from the Gaussian source term in the London 
equatio: 



6,17 



p(r) 



^0 



exp 



(8) 



Although this form is presumed to be valid only in the GL 
limit (i.e., near Tc), it is known that Eq. ([7j) reproduces 
the cutoff function derived from the most commonly as- 
sumed form of the order parameter that fits well the so- 
lution of the Bogoliubov-de Gennes (BdG) equation. 



V^oo tanh-. 



(9) 



We also note that a renormalzation factor, Vl — b (where 
b = B/Bc2)^ was introduced to the m-London model in 
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FIG. 1: Density distribution function {n{B), left) and cor- 
responding spatial field profile (B(r), right) for a hexagonal 
flux line lattice obtained from the Ginzburg-Landau model; 
A = 125 nm, ^ = 30 nm, and Bo = 0.15 T. The distance 
between flux lines, ao, is 12.6 nm. The points marked as 
'C', 'S', 'V in n{B) respectively correspond to those in B(r). 
Thus, despite its one-dimensionality, n{B) has a site-selective 
sensitivity in the unit cell of FLL. 




the earlier literatures to account for the reduction of the 
GL order parameter at higher fields (see below), where A 
and in Eq. (|6|) are divided by this factor. However, 
it turns out throughout this work that this renoramliza- 
tion is too strong to be consistent with our experimental 
observation, and therefore disregarded in the following 
analysis. 

On the other hand, an approach using variational 
method has been made to obtain the analytical GL solu- 
tion for B{r) by assuming 



which leads to a squared Lorentzian source teim^ 



Pir) 



(10) 



(11) 



and associated cutoff function (for B <C Bc2) 



F{K,^,) = uK^{u), u = ^/K^^T^J}^, (12) 

where Ki{x) is the modified Bessel function and 
V^^o is the variational parameter This model was ex- 
tended to higher fields and to anisotropic order param- 
eters by Hao et al}^ to consider the reduction of order 
parameter due to the overlap of vortex cores, yielding the 
Fourier component 



5(K)=i.o^, u 

uKi{v) 



VK^^[+^,v=^foo, (13) 



where and foo are the variational parameters with the 
latter representing the field-dependent order parameter 



FIG. 2: Superconducting transition in YBe observed by (a) 
magnetic susceptibility {/j^qH = 1 mT), and (b) the upper 
critical field determined by specific heat. Solid curve in (b) is 
fit by a power law (see text). 



(V^oo /ooV^oo; foo 1 for the dilute vortex limit). Un- 
der a condition of A^i^^-^^ 1, where Kmin is the small- 
est nonzero reciplocal lattice vector, Eq. (p!3|) is simplified 
for the extreme type II case 1) using an approxi- 

mation Ki{x) 1/x to yieldi^ 



b{K) 



2 uKiju) 



where 



b = B/Bc2 



26(1 - bf 



(14) 



(15) 
(16) 



In the present case of YBe, ^min for a hexagonal FLL 
is 4>/3/ao — 0.015 nm~^ (where ao — 455 nm for the 
lowest field Bd ^ 0.02 T), which leads to A^iT^^^ 3.9 
(A :^ 130 nm), indicating that A^iiT^-^ is considerably 
greater than unity. Moreover, i^^in^^ — ^-^ — nm), 
which is again larger than v'^ < 0.05 to justify 

the approximation, u K(_y [see Eq. (p!3|) ]. Thus, YBe 
satisfies the conditions for Eqs. (p!4|) -(p!6 |) to be employed 
as a model to describe B{t). 

An example of B{t) and corresponding n{B) calcu- 
lated by the GL model is shown in Fig. [TJ Regardless of 
the details in the modeling of 5(r), it predicts an asym- 
metric field profile for n{B) characterized by a negatively 
shifted sharp peak due to the van Hove singularity asso- 
ciated with the saddle points of B{r) (marked as 'S' in 
Fig.[T]), and an adiabatic tail towards higher fields where 
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FIG. 3: Real amplitude of the fast Fourier transform (FFT) of 
/iSR time spectra under several temperature/magnetic fields 
(/jqH) in the mixed state of YBe (corresponds to n(B)). The 
small bump at B — /j^qH = seen in the spectra at low temper- 
ature/field is due to background signal from muons stopped 
in the sample holder. 



the maximal field is determined by 5(|r — r^l ~ The 
nearly one-to-one correspondence between the character- 
istic points of n{B) and those on B{r) makes it feasible 
to deduce complex physical parameters including A and 
^ from a single /iSR spectrum. Meanwhile, information 
on the structure of FLL must be provided by other tech- 
niques (e.g., small angle neutron scattering). Since no 
such information seems to be available for YBe, we pre- 
sume in the following that a hexagonal FLL is realized. 

While one might argue that the m-London model 
should be replaced by the GL model simply because of 
the wider field range of applicability, we point out that 
the physical meaning of in the GL model has some 
uncertainty in the sense that it has a complicated field 
dependence as a cutoff parameter [see Eq. (p!6]) ]. yield- 
ing different values at different fields for a common ^. 
This casts some ambiguity, particularly on the behavior 
of vortex core radius as a function of external field. The 
GL model has another disadvantage in the practical ap- 
plication that it requires significant computing resources 
(CPU power, in particular) for the precise treatment. 
The m-London model, in contrast, is relatively free from 
such ambiguity of model parameters, and it can be im- 
plemented with ordinary computing environment. Thus, 
the m-London model would remain to be a useful basis 
for the application of /iSR to the study of the mixed state, 
as far as precaution is taken on the limit of its validity. 

Finally, it may be worthy of mention that F(i^, ^^) 
would be also dependent on the details of the Fermi sur- 
face, and thereby it may be subject to further correc- 
tion to include K dependence.^^ While the effect of such 
anisotropy within the vortex cores may be small at the 
dilute vortex limit, it may be important at fields close to 
Bc2 in anisotropic superconductors. 



III. EXPERIMENTAL DETAILS 

Single crystals of YBe were prepared by the floating 
zone technique, where the details of the crystal growth 
is found elsewhere. Magnetization, specific heat, and 
electronic resistivity were measured prior to /iSR experi- 
ment to evaluate the bulk properties of the specimen. As 
shown in Fig.O the result indicates that the crystal has a 
superconducting transition temperature of Tc = 6.95(15) 
K (defined at the midpoint of change in xo) and an upper 
critical field of 5c2(0) = 0.31 T. Other physical proper- 
ties deduced from those measurements are summarized 
in Table H The ratio 2A(0)/A:bTc 3.67 (deduced from 
the jump of the electronic specific heat at Tc) is close to 
the BCS value of 3.53. The mean-free path is estimated 
to be ~2.83 nm from the residual resistivity ratio. This is 
much shorter than [= (^o/27r5c2(0))"^/^], thereby in- 
dicating that the superconductivity is in the dirty limit. 
This strongly suggests that nonlocal effect in terms of K 
dependence, if it exists at all, must be masked by the 
impurity scattering in this specimen. 

Conventional /iSR measurements under a transverse 
field (TF) were carried out on the M15 beamline at 
the Tri-University Meson Facility (TRIUMF, Vancouver, 
Canada), where the field direction was perpendicular to 
the crystal plane (normal to [001] axis) covering an area 
of 8 mm X 8 mm on the sample holder {H \\ z and [001]). 
A positive muon beam with the initial polarization nor- 
mal to the z axis was implanted into the specimen, and 
the time evolution of muon spin polarization [P{t)] was 
monitored by detecting the decay positrons emitted pref- 
erentially towards the direction of P{t). Each measure- 
ment was done by accumulating 2 ~ 4 x 10^ positron 
events on two pairs of scintillation counters (for x and y 
directions) that yield /iSR signals 

A{t) = A,{t) + iAy{t) = AoP{t), (17) 

where Ax{t) and Ay{t) are the decay positron asymme- 
try for the respective pair of counters with Aq being the 
total asymmetry. Measurements on the vortex state were 
made under field-cooled conditions to minimize the effect 
of flux pinning that may give rise to additional broaden- 
ing of the /iSR lineshape. 



Tc (K) rrr B^i (mT) Bc2 (T) / (nm) Al (nm) (nm) 

6.95(15) 3.67 4.62 18.5 0.31 2.83 134(2) 33(1) 



TABLE L Bulk properties of YBe used for the /iSR experi- 
ment, where Tc is the superconducting transition temperature 
(defined by the midpoint of xo), A(0) is the energy gap ex- 
trapolated to zero temperature, rrr is the residual resistivity 
ratio, Bel and Bc2 are the lower and upper critical field, / 
is the mean free path (deduced from rrr), Al is the London 
penetration depth (from Bd), and is the BCS coherence 
length (from Bc2)- 
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FIG. 4: Time revolution of /x-e decay asymmetry [Ax{t)^ cir- 
cles] in YBe at iiqH = 0.15 T and T=2.0 K displayed on 
a rotating reference frame of 18 MHz, where solid curves in 
a) and c) are fits assuming B{y) obtained by the m-London 
model (A — 153 nm, — 25.6 nm) and GL model (A — 123 
nm, ^ — 29.6 nm), respectively. Fit errors are shown in b) 
and d) as difference patterns between data and fit. 



IV. RESULT 

A. Comparison of models: Field dependence of A 
and ^ 

Examples of the fast Fourier transform (EFT) of /iSR 
spectra at several external magnetic fields are shown in 
Fig. O where one can readily observe the asymemtric 
lineshape with slight smearing due to the limited time 
window for EFT (0-6 /is). The peaks at lower field rep- 
resent the van Hove singularity due to the saddle points 
of B{v) ('S' in Fig. [ij). additional peak at 5 - iiqH = 
is due to muons which missed the specimen and stopped 
in the positron counters surrounding the specimen. The 
value of jiioH is experimentally determined by the spec- 
trum of each field above Tc. Based on the least-square 
method with appropriate consideration of the statistical 
uncertainty, the /iSR spectra are compared (in time do- 
main) with those calculated by using Eqs.([2])-(|6|), 

/oo 
[(1 - /b)e-K+-n)*^n(i?) 
-OO 

+ fkSiB - iioH)y^^''*-'UB, (18) 

from which a set of parameters for n{B)^ A, and [or 
^ in Eq. ([16]) for the GL model] is deduced. Here, <jp is 
the additional relaxation due to random flux pinning, 
is that due to random local fields from nuclear magnetic 
moments, and /b denotes the fractional yield of the back- 
ground signal. (Jn (— 0.17 MHz) was determined from the 
/iSR spectrum above Tc, and /b was less than 3 % of the 
total asymmetry. The background comes from muons 
stopping in the positron counters (which survived due to 
the finite deficiency of veto logic circuits to eliminate such 
positron events) and thereby showing negligibly small re- 
laxation rate. We found that the average internal field 
(Bq) showed a small positive shift (< 1 mT) from the 
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FIG. 5: Residual x divided by the number of degree of free- 
dom (A^df)- 



external field {ijlqH) with randomly varying magnitude. 
We attribute this to the effect of random pinning of vor- 
tices and associated demagnetization which seems to be 
strong in this specimen as suggested by relatively large <jp 
(!::^0.2-0.5 MHz at any field/temperature). Thus, we al- 
lowed ^0 to vary slightly (< 1 mT) to obtain the best fit. 
Fig. m shows an example of analysis in the time domain, 
where fits by the m-London and GL models are displayed 
for the data at iiqH = 0.15 T and T = 2.0 K. Although 
both models yield reasonable fit to data, the ratio of 
to the number of degree of freedom is considerably re- 
duced from 1.55 for the m-London model to 1.39 for the 
GL model: one may observe the scatter in b) is slightly 
reduced in d). As is evident in Fig. [5l this tendency is 
more clearly observed at higher fields, and it is consistent 
with the presumption that the GL model provides bet- 
ter description over the higher field region. Meanwhile, 
it should be also stressed that both models yield almost 
the same quality of fits at lower fields {b < 0.3). 

Figure [6] shows comparison between those two models 
for the profiles of 5(r), supercurrent calculated by the 
Maxwell's relation. 



J(r) = cm\{B{r)} = 



dB,{r) dB,{r) 



dx 



dy 



(19) 



and superconducting order parameter [?/^(r)/?/^(oo)] plot- 
ted along the straight line connecting the nearest- 
neighbor vortices, where the curves in (a)-(b) and (d) 
are calculated by using parameters obtained by fits shown 
in Fig. m The field profile generated by the GL model 
keeps increasing towards the vortex center beyond that 
by m-London model, yielding smaller effective radius (ro) 
for vortex cores defined by the peak of the supercurrent 
[^(^o) = 1^1 max with ro — 25 nm for GL and 30 nm 
for m-London]. Note that the definition of ^ in the GL 
model is quite different from in the m-London model, 
and thus they cannot be compared directly. Meanwhile, 
ro is defined directly from the field profile using Eq. (p!9|) , 
and thereby it serves as a common ground for compar- 
ison of vortex core size. It is noteworthy that, despite 



6 




20 40 60 80 100 120 

(nm) 



-2 - 1 



FIG. 6: (a) The spatial field profile of B{r) along the straight 
line connecting nearest neighboring vortices [shown in (c) 
by dashed line], where A and (0 taken respectively 
from the result of fits shown in Fig 2] (b) The spatial 
profile of the supercurrent obtained by the Maxwell's rela- 
tion [J{r) = curlB(r)]. The normalized order parameter 
[ip(r) = ip{r)/ip{oo)] calculated by Eq. (|lQp is also plotted, (d) 
The density distribution n{B) corresponding to each curve in 
(a). 



relatively small value of b 0.5), the order parameter 
in Fig. [6jb) calculated as a product of Eq. (p!Q|) for two 
neighboring vortices exhibits considerable reduction at 
the center between two vortices. This again suggests that 
the GL model must be used for the proper description of 
FLL state over the relevant field range. 

While the difference between two models appears to 
be localized in the field profile near the vortex cores, it 
gives rise to change in the deduced value of A. As shown 
in Fig. [71(a), those obtained by the m-London model ex- 
hibits a clear tendency of longer A at higher fields than 
m-London model. A similar result is reported for the 
case of NbSe2 (see Fig. 32 of Ref. Q), although the data 
are limited to lower fields {b < 0.31). In the case of YBe, 
this leads to qualitative difference in the dimensionless 
paramater, 77, to describe the gradient of A against field, 
which is defined as 
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FIG. 7: Magnetic field dependence of penetration depth and 
cutoff parameter as determined from the m-London and the 
GL models. The solid and dashed lines are linear fits to the 
data obtained by the respective models. 



ing the magnitude of r] to employ model appropriate for 
the relevant field range of data. We will come back to 
this issue later below. 

In contrast to the case of A, there is not much differ- 
ence in the field dependence of cutoff parameters between 
two models, where both and ^ fall into similar values 
between 20 and 30 nm. While ^ shows a weak tendency 
of decrease with increasing field, does not show any 
clear trend. As mentioned earlier, the corresponding vor- 
tex core radius defined as the peak position of supercur- 
rent is slightly smaller than ^ in the GL model, while an 
opposite tendency is observed for the case of m-London 
model. Taking an average of ^ over the observed field 
range [= 27(2) nm], we obtain the corresponding GL pa- 
rameter n = 4.7(4) in good agreement with n A.l from 
the bulk property measurements. The fact that both A 
and ^ are independent of H means that n is also a con- 
stant against H. 



A(6) = A(0)[l + 7?-6]. (20) 

Fits of data in Fig. [T^a) using the above relation yields 
T] = 0.59(3) for the m-London model and rj = —0.01(3) 
for the GL model. Since it is anticipated that 7^ for 
the case of isotropic order parameter, the analysis by the 
GL model, which provides better description of 5(r), im- 
plies that YBe belongs to the class of conventional BCS 
superconductors with isotropic energy gap. This in turn 
indicates that a precaution must be taken upon evaluat- 



B. Temperature dependence of A and ^ 

Figure [8] shows the temperature dependence of 1/ A^ 
and ^ determined by the fits of /iSR spectra at (IqH = 
0.05 T, using solutions of the GL model for generating 
n{B). Since we showed in the previous section that the 
GL model is superior to m-London model, we focus on 
the result of analysis by the GL model in the following. 

The solid curve in Fig. [81(a) is a fit by the weak coupling 
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BCS theory 
1 1 



A2(T) A2(0) 



\ EdE 
dEJ _ A2 



, (21) 



where f{E) is the Fermi distribution function and A = 
A(T) is the temperature-dependent gap energy. It re- 
produces data reasonably well, yielding A(0) = 123(2) 
nm and Tc = 6.2(1) K. The obtained Tc is in excel- 
lent agreement with that determined by bulk measure- 
ments at iiqH = 0.05 T shown in Fig. [2fb). The scat- 
ter of data points might be attributed to the relatively 
strong pinning that might give rise to a distortion of B{y) 
that is beyond the limit of Gaussian approximation [i.e., 
exp(— cTpt^), see the inset of Fig. [HKa)]. Although the 
absence of data below 2 K does not allow definite conclu- 
sion, the temperature dependence of 1/A^ is consistent 
with the conventional BCS model for isotropic 5-wave 
pairing, supporting the conclusion drawn from the field 
dependence of A. 

The temperature dependence of ^ is displayed in 
Fig. mb). It exhibits a tendency of decrease with de- 
creasing temperature and saturates below T < 4 K. Sim- 
ilar results are reported for the cases of NbSe2 (Refl2ll) 
and V (RefilO), where the observed temperature depen- 
dence has been attributed to the so-called Kramer-Pesch 
effect. In the FLL state, the order parameter serves as 
a potential energy [ex V^(r), the pair potential] for the 
quasiparticles, and they are trapped in the vortex cores 
(by the Andreev refiection) to form discrete energy levels. 
In clean superconductors, it is then predicted that the 
vortex core radius (ex ^) must be dependent on tempera- 
ture, because the higher energy levels, which are occupied 
by quasiparticles at higher temperatures, have greater 
extension corresponding to a larger core size. Unfortu- 
nately, however, the experimentally observed tempera- 
ture dependence of ^ seems to be much weaker than the 
predicted linear dependence, ^ ~ £^o-{T/Tc)^ in those ear- 
lier examples which are in the clean limit. Meanwhile, 
it is predicted that the effect is sensitive to the impu- 
rity scattering and thereby it would be smeared out in 
the dirty limit corresponding to the present case of YBe- 
This is again in contrast to the observed behavior of ^ 
for T > 4 K, suggesting that it cannot be attributed to 
the Kramer-Pesch effect. 

Here, we recall a much simpler argument that the 
weak temperature dependence of ^ can be naturally ex- 
plained by considering the fact that Bc2 varies with 
temperature;^ Since the GL coherence length near Tc 
is given by the relation 



eGL(T) 



^0 



27rBc2{T)' 



(22) 



it is natural to assume ^(T) :^ Cgl(T) oc [5c2(T)]-^/^ 
Thus, the increase of ^ at higher temperatures might be 
attributed to the decrease of Bc2{T). A fit of Bc2{T) in 
Fig. [2](b) by a power law. 
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FIG. 8: Temperature dependence of 1/A^ (a) and cutoff pa- 
rameter (b) as determined from the GL model. Solid curve in 
(a) is a fit by the weak coupling BCS theory, and that in (b) 



is fit by a power law proportional to [Bc2{T)] 



-1/2 



text). 



B,2{T) = B,2mi-aT/Tr], 



(23) 



Inset: The rate of additional relaxation due to flux pinning. 



yields Bc2{0) = 0.314(1) T, C = 0.970(1) and u = 1.56(1) 
The solid curve in Fig. ^h) is the best fit obtained by 
assuming Eq. ([22|) and ([23|) with the above parameter 
values, which yields Tc = 7.6(4) K. The agreement is 
satisfactory and thereby it demonstrates that ^ at 2 K 
is determined by Bc2 near Tc. This indicates that the 
thermal fiuctuation is relatively strong over the relevant 
temperature region (T/Tc > 0.29), and that the Kramer- 
Pesch effect must be searched for at temperatures low 
enough so that one may be able to neglect the tempera- 
ture dependence of Bc2{T). 



V. DISCUSSION 

A. Superconductivity of YBe 

Despite the early discovery of superconductivity in 
YBe, there are not much literatures published so far on 
its superconducting property. Fortunately, a recent pa- 
per by Lortz et al. reports detailed measurements on spe- 
cific heat, resistivity, and thermal expansion.^ The spe- 
cific heat in the sperconducting state is that typically 
found for a single-band, isotropic BCS superconductor, 
which we also confirmed on our specimen by similar mea- 
surements. Meanwhile, the electron-phonon interaction 
turns out to be much stronger than that in other boride 
superconductors such as ZrBi2 and MgB2. They report 
an enhanced coupling to phonon modes lying near 8 meV 
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in YBfi, which is relatively close to ~ 15 meV of ZrBi2 
(Ref i24) . This is to be compared with ~ 60 raeV in 
MgB2 (Ref I25I) . thus it provides an explanation for the 
difference of Tc among those three borides. Interestingly, 
those low-energy frequency modes are attributed to the 
vibration of Y or Zr atoms loosely bound in oversized 
boron cages that is now called "rattling" motion. Here, 
the longer metal to boron bond length in YBe leads to 
a weaker force constant and larger vibrational amplitude 
favorable for higher Tg. The reduction of Tc under high 
pressure predicted from the thermal expansion measure- 
ment, which is also anticipated from the coupling to the 
rattling phonon, has been confirmed by the recent exper- 
iment on the pressure effect 

In the previous section, we have shown that the field 
and temperature dependence of A is fully consistent with 
the presence of isotropic order parameter in YBg. It 
does not show any hint of nonlocal effect nor of low-lying 
quasiparticle excitation associated with the anomaly in 
the order parameter. This in turn suggests that the 
strong coupling of electrons to the rattling phonons does 
not necessarily lead to anomaly in the superconducting 
order parameter. The situation is readily understood 
by considering that YBe has a highly symmetric three- 
dimensional (3D) band structure, while MgB2 is char- 
acterized by two-dimensional a band and 3D-7r band 
which lead to multi-gapped superconductivity (see be- 
low). However, it should be noted that the present spec- 
imen is in the dirty limit where the anisotropic feature 
having the length scale longer than the mean-free path 
(/ = 2.83 nm, corresponding to Hvy/ttI ~ 7 meV assum- 
ing that ~ 10^ m/s) must be smeared out by scat- 
tering of electrons. Therefore, further study with much 
better crystal quality may be needed for the definitive 
conclusion. 



?7 as a criterion of anomaly in the order 
parameter 



We saw in the section HV Al that the slope 77 of the field 
dependence of A obtained by m-London model may be 
slightly overestimated when the model is applied beyond 
the presumed limit, B / Bc2 ^ 0.25. This may bring about 
a concern on the past literatures discussing the degree of 
gap anisotropy by the magnitude oi rj. To make the sit- 
uation clear, we collected results reporting the value of r] 
obtained by using m-London model, which are shown in 
Fig. [9] as a plot of 77 versus the maximal field of measure- 
ment (/io^max) normalized by fioHc2 at the measured 
temperature. As a comparison, those obtained by the 
GL model is also plotted for VaSi (Ref. 27) and YBe (the 
present result). 

It is relatively well established by other experimen- 
tal techniques that the group of compounds shown by 
the filled symbols in Fig. [9] (those with r] > 1) have 
anisotropic superconducting gap or multi-gapped, ex- 
cept for the case of KOsaOe (RefsHHl) in which the 
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FIG. 9: Slope of field dependent A plotted against the max- 
imal field of measurement normalized by the upper critical 
field, where the m-London model was used to generate B{r) 
(the relaxation was approximated by a Gaussian damping for 
CdsOssOr and MgBs) [(Refs|8iL23:27:28:29:30:3L 32^333- 
Filled symbols: superconductors with anisotropic gap (cir- 
cles) or multigapped (squares). Open circles and triangles: 
those with isotropic gap (where the GL model was used for 
the latter case). 



situation is not clear at this stage. While the obser- 
vation that 77 > 1 for those compounds supports the 
use of 7^ as a criterion for the gap anisotropy, the cases 
of YNisBsC (ReflU) and MgBs (Refill) might require 
closer examination because of large i^max/^c2- Fortu- 
nately for the case of YNi2B2C, data with small field 
step exist at lower fields (h < 0.2), and one can clearly 
observe that 77 remains to be 1 within the limit of m- 
London model {h ^ 0.25). This may be partly attributed 
to the presence of point nodes in the superconducting 
order parameter^^^^^ that would give rise to the nonlo- 
cal effect in YNi2B2C at higher fields. Since the nonlocal 
effect directly affects the behavior of A, it may be less sen- 
sitive to the difference in the modeling of vortex cores be- 
tween the m-London model and GL model. In particular, 
the situation might be that the magnitude of rj at higher 
fields, which is predominantly determined by the nonlo- 
cal effect, coincides with that due to the Volovik effect 
observed at lower fields. Interestingly in YBa2Cu306.95 
(YBCO, having line nodes), it is reported that 77 2 for 
h > 0.04 which is much smaller than 77 > 5 for 6 < 0.02 
(Reflsl), where the magnitude of 77 at higher fields is ex- 
plained by considering the nonlocal effect associated with 
the line nodes. A similar situation would be anticipated 
with less magnitude of r] at higher fields for the case of 
YNi2B2C with point nodes. The point node also provides 
a natural explanation for small rj compared with that of 
YBCO at lower fields. 

Here, it might be stressed again that there are two dif- 
ferent sources of nonlocal effect, the one associated with 
the anisotropy in the superconducting order parameters 
and that associated with the anisotropy in the Fermi sur- 
face. In general, the latter effect manifests itself as the 
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transformation of FLL from hexagonal to squared with 
increasing external field, where transformation occurs at 
relatively low fields (e.g., b < 0.06 in YNi2B2C). Mean- 
while, the anisotropy in the order parameter has rela- 
tively smaller energy scales and thereby it emerges at 
higher fields. While both yield a qualitatively similar 
niodification to the effective London penetration depth 
(L), their consequences are independent (even compet- 
ing) with each other 

It is now well established that MgB2 is characterized 
by two different energy gaps corresponding to its band 
structure (A^r — 6.8 meV for a band and A^^ 2.2 meV 
for TT band).^^ The analysis in Ref.29 was made before the 
revelation of the double gap, using the Gaussian approx- 
imation instead of m-London model. However, it cap- 
tures anomalous behavior of (proportional to the 
Gaussian linewidth) unexpected for the conventional su- 
perconductors with single gap. 

The Gaussian model is based on a crude approximation 
that the Gaussian damping rate, a, is primarily deter- 
mined by the penetration depth (valid for extreme type 
II superconductors), 

a = 7m(EW)'^'cxG(5)A-2, (24) 

K 

G(6) ^(1-6)[1 + 3.9(1 -6)Y/', 

where the function G{b) represents the reduction of a 
mainly due to the overlap of vortices (proportional to 
1 — b) and additional narrowing due to the contribution 
of vortex cores. ^ The presence of two gaps means that 
there are two different upper critical fields corresponding 
to respective coherence length, 

^c2(i) = ^;^, {i = cr,7T), (25) 

from which one may anticipate that 

a (X [a^G{b^) + (1 - a,)G(6,)]A-^ (26) 

where bi = /io^/^c2(i) (with 5c2(7r) < ^c2(a)), «a is the 
fractional weight of a component, and G{bT^) must be set 
to zero for b^^ > 1. The second term in the above equa- 
tion, which exhibits steeper field dependence than the 
first term, strongly reduces a with increasing field over 
the field range of < jiioH < 5c2(7r)- From the viewpoint 
of single Bc2 , such behavior appears as an anomalous de- 
viation of A from Eq. ([24|) and thereby easily identified as 
an increase of 77 as reported in Refl29l. As a matter of 
fact, recent attempt to analyze a assuming a model sim- 
ilar to Eq. ([261) reports = 5.1(2) nm and = 23(1) 
nm.— Considering that the magnitude of is as large 
as nearly a half of A [=49(4) nm], the enhancement of 
A is physically understood as that due to the QP's over 
a radial region < r < .^7^ around vortex cores. This 
is presumed to be a common situation for multi-gapped 



superconductors including another example of NbSe2 4£ 
We also point out that the smaller ^(T) (e.g., ^7r(T) in 
MgB2) may exceed A at a certain temperature, as it in- 
creases with T [see Eq. (|22]) and Fig [81(b)]. Then, there is 
a possibility that the nonlocal effect may set in to further 
enhance the effective A. 

Now it is clear that the magnitude of r] indeed serves 
as a criterion for the anisotropy of superconducting order 
parameter, where that obtained from the analysis using 
the m-London model remains useful (particularly those 
for which the maximal field range is within the presumed 
limit of validity, b < 0.25). In this context, the previ- 
ous result for KOS2O6 may be regarded as a strong case 
for the presence of anisotropic order parameter (includ- 
ing multi-gap) in this compound. Considering results of 
other experimental techniques that the order parameter 
is nodeless^ further analysis of /iSR data by a model 
based on the double-gap scenario is now in progress. 



VI. SUMMARY AND CONCLUSION 

We showed that the effective London penetration 
depth in YBe, obtained by analyzing /iSR data using 
the field profile generated by the analytical GL model, 
is independent of external field {r] c:^ 0) over a region 
< iiqH/Bc2 < 0.65. We made a comparative analysis 
using the modified London model and GL model, from 
which it was inferred that the GL model should be used 
for the precise determination of 77 when the field range 
of measurements extends over the presumed limit of va- 
lidity for the m-London model {iiqH/Bc2 > 0.25). We 
also showed that the temperature dependence of 1/A^ is 
perfectly in line with that predicted for isotropic BCS su- 
perconductors. This, together with the observation that 
77 0, supports for the conclusion drawn from bulk mea- 
surements that YBe belongs to the class of conventional 
BCS superconductors with isotropic s-wave pairing. Fol- 
lowing this result, we re-exmined the past reports on r] 
obtained by the analysis based on the m-London model 
(and those based on the Gaussian approximation), and 
confirmed that the conclusions drawn for respective com- 
pounds remain unchanged. Thus, the m-London model 
continues to be a useful ground for the /iSR study of 
vortex state within the limit of its validity. 
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